Bounds on fifth forces from precision measurements on molecules 
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Highly accurate results from frequency measurements on neutral hydrogen molecules H2, HD 
and D2 as well as the HD^ ion can be interpreted in terms of constraints on possible fifth-force 
interactions. Where the hydrogen atom is a probe for yet unknown lepton-hadron interactions, and 
the helium atom is sensitive for lepton-lepton interactions, molecules open the domain to search 
for additional long-range hadron-hadron forces. First principles calculations in the framework of 
quantum electrodynamics have now advanced to the level that hydrogen molecules and hydrogen 
molecular ions have become calculable systems, making them a search-ground for fifth forces. Fol- 
lowing a phenomenological treatment of unknown hadron-hadron interactions written in terms of a 
Yukawa potential of the form V^i^r) — /3exp(— r/A)/r current precision measurements on hydrogenic 
molecules yield a constraint /? < 1 x 10~^ eV-A for long-range hadron-hadron interactions at typical 
force ranges commensurate with separations of a chemical bond, i.e. A « 1 A and beyond. 

PACS numbers: 



I. INTRODUCTION 

While the Standard Model (SM) of physics explains 
physical phenomena observed at the microscopic scale, 
phenomena of Dark Matter [1] and Dark Energy [2] at the 
cosmological scale are considered as unsolved problems, 
possibly hinting at physics beyond the SM. String the- 
ory [3] and supersymmetry [4] seek to accomodate these 
phenomena, as well as gravity, with the SM within a 
unified model. New kinds of fundamental interactions 
and/or extra dimensions are postulated as extensions of 
the SM |5l IS , which could be probed via high-energy par- 
ticle colliders. However, there is also a frontier of low- 
energy physics [7j , with predictions of weakly interacting 
particles at the eV energy scale, that could be probed in 
table-top experiments. 

Celebrated examples of probing new physics at the 
atomic scale include experiments aimed at measuring 
an electric dipole moment of the electron, first through 
measurements on atoms [H], and later at increased pre- 
cision level also on molecules [5]. Experimental searches 
have been carried out to detect anomalous spin-spin in- 
teractions between electrons investigated through spec- 
troscopy in ion traps [T^, or in paramagnetic salts [TT] . 
Similar anomalous spin-spin couplings between neutrons 
are investigated in '^He/^^^Xe masers [12] . Constraints on 
short-range spin-dependent interactions between protons 
and deuterons at the A scale are derived from nuclear 
magnetic resonance experiments in the HD molecule (15] . 

To date the energy level structure of atomic and molec- 
ular systems can be fully described by electromagnetism, 
in its most advanced form by the theory of quantum elec- 
trodynamics (QED). Effects of the weak force, leading to 
parity non-conservation, have been clearly observed at 
the atomic scale [Hj, but not yet in light atoms as H/D 



and He. In these calculable systems effects of the weak 
interaction on the energy level structure is orders of mag- 
nitude away from experimental determination, although 
for muonic hydrogen the contribution is just below the 
accuracy of the proton-size contribution [15] . Strong in- 
teractions, in quantum chromodynamics confined to the 
fm scale, are at the origin of nuclear g^ factors and nu- 
clear spins In, thus influencing atomic and molecular 
level energies in terms of hyperfine structure. Obviously, 
gravitational interactions are far too weak to play any 
role in the calculation of level energies in atomic sys- 
tems. Hence the rationale for probing fifth forces beyond 
the SM is based on a search for deviations from QED in 
the quantum level structure of calculable systems at the 
atomic scale. 

QED has been tested in light atoms to high precision, 
for example by the observed agreement in the derived 
proton rms charge radius rp from several atomic hydro- 
gen transitions. In fact, by assuming the correctness of 
QED, the present CODATA recommendation fW for Tp 
is based on these derived values from high-precision hy- 
drogen spectroscopy along with the value from electron- 
proton scattering experiments, where both derivations 
are in good agreement. However, this rp value is in 
7a-disagreement with that derived from recent measure- 
ments on muonic hydrogen (^~p+) [TT] [18]. Although 
this deviation is interpreted as a puzzle on the proton 
size, the solution to this conundrum might as well be 
found in unaccounted effects within QED ,19i . 

The He atom, a two-electron system, is still accessible 
for accurate ab initio calculations of the level structure in 
the framework of QED [2D1 [H] , although measured level 
splittings arc in several cases more accurate than the- 
ory [22H24] . Precision measurements on the helium atom 
allow for tests of QED including lepton-lepton interac- 



tions that are not detectable in measurements of e"!!"*", 
nor in anti- hydrogen (e"*"!!"), nor in muonic hydrogen. 

Precision QED test have been extended to molecules, 
i.e. to systems with more than one nucleon, where 
long-range hadron-hadron forces come into play. Quan- 
tum ab initio calculations have been performed on the 
single-electron Hj ionic system by solving the three- 
body Coulomb problem and calculating QED contribu- 
tions up to high order in a [23 [2S] ■ Similar calculations 
of equal accuracy have been carried out for the HD+ 
ion [27H5S], which exhibits a small dipole moment and 
is therefore more amenable for laser precision measure- 
ments [SniEI]. In recent years, full-fledged level structure 
calculations, including QED and high-order relativistic 
contributions have also been carried out for the neutral 
hydrogen molecules and the deuterated isotopomers [321" 
[34] . These calculations have been tested in the deter- 
mination of dissociation limits of H2 "36], D2 [37] and 
HD 38J , as well as in measurements of the ground state 
rotational sequence of H2 |39| , and the vibrational split- 
tings in H2, D2 and HD [40H44] . 

Fifth-force tests are commonly associated with testing 
non-Newtonian gravity, sometimes motivated by theories 
postulating extra dimensions [S]. Tests of the inverse- 
square law behaviour of gravity have been carried out 
over a wide distance scale from kilometers to submicron, 
where the present short distance constraints are obtained 
from Casimir force experiments |45j . We propose here 
that molecules open up a new arena for probing fifth 
forces at typical distance separations occurring in chem- 
ical bonds, thus at length scales naturally set to A dis- 
tances. Attractive or repulsive additional forces can be 
probed from precision metrology measurements on cal- 
culable molecular systems. Based on recent precision 
measurements on HD"*" ions and H2, D2 and HD neutral 
molecules, and comparison with advanced QED calcula- 
tions, constraints on a fifth-force can be derived, which 
is parameterized by a generalized Yukawa potential for 
a certain effective range. Where lepton-nucleon interac- 
tions and lepton-lepton interactions may be probed in 
atomic hydrogen and helium, in the present study addi- 
tional forces between hadrons at long range are targeted, 
for which molecules are a good test ground. 



II. QED CALCULATIONS IN NEUTRAL AND 
IONIC MOLECULAR HYDROGEN 

Accurate ab initio level energies E of molecular hy- 
drogen and its deuterated isotopomers are calculated in 
the framework of nonrelativistic quantum electrodynam- 
ics (NRQED) from an evaluation in orders of the electro- 
magnetic coupling constant a 



E{a) = f (") + a" £:(2) -f a^ gi^) + a^ f W 



(1) 



The nonrelativistic energy £^°'' is obtained by solving the 
Schrodinger equation using the Born-Oppenheimer po- 
tential with 15 digit accuracy [46] . Adiabatic and nonadi- 



abatic corrections are subsequently calculated perturba- 
tively in powers of the electron-nucleus mass ratio. This 
procedure results in nonrelativistic binding energies ac- 
curate to a few parts in 10""' cm"-'^ [37]. The u — g 
symmetry-breaking is taken into account for the specific 
case of HD [34] . Leading-order relativistic corrections S^^'^ 
have been calculated from the expectation value of the 
Pauli Hamiltonian 32 , 33 . The leading QED corrections 
of order f ^^^ are treated similarly as for the hydrogen and 
helium atoms [48]. This results in an accuracy of 10~^ 
cm~^ for £^^^ in molecular hydrogen due to neglect of 
nonadiabatic and relativistic recoil corrections. The main 
contribution to the calculation uncertainty at present is 
from the f ^^^ QED correction, which has not been calcu- 
lated explicitly due to the high complexity of NRQED op- 
erators. The radial nuclear functions for v = and v — 1 
probe almost the same range of internuclear distance, 
leading to significant cancellation of the uncertainty of 
the energy contributions. The final theoretical predic- 
tions are estimated to be accurate to 1 x 10"* cm"^ for 
the full QED evaluation of the rotationless fundamental 
ground tones in H2, D2 and HD. There is less cancella- 
tion in uncertainty of transition energies as the difference 
in the vibrational quantum number Av increases, lead- 
ing to larger uncertainties at ^ 1 x 10"^ cm"^ for the 
dissociation energies {Dq) of molecular hydrogen and its 
isotopomers. 

A reduction in complexity of three-body system molec- 
ular ions compared to the (two-electron) neutral molec- 
ular species results in better calculation accuracy for Hj 
and HD+. Nonrelativistic energies are obtained with up 
to 30-digit numerical precision for the low-lying vibra- 
tional states [15] using a direct variational approach. Rel- 
ativistic and QED corrections are also calculated using 
the NRQED framework similarly expressed as Eq. (fl]) for 
the neutral species. The leading-order relativistic cor- 
rections f'^-* have been calculated with sub-kHz accu- 
racy |50j . A recent improvement in the evaluation of the 
Bethe logarithm [2S] has enabled an increased accuracy 
of the 5*^^^ term at better than 50 Hz. The total uncer- 
tainty is dominated by the contribution of the higher or- 
der QED terms £^'^^ and £'^^\ estimated to be ~ 20 kHz. 
Note that in the case of the ion with its unpaired elec- 
tron, the hyperfine interaction (which is absent to first 
order in the neutrals) is addressed separately [51]. The 
QED calculation uncertainty contributes 21 kHz while 
the (in)accuracy of fundamental constants contributes 10 
kHz to the total uncertainty of the ab initio calculation 
for the HD+ w = — >■ 1 R(l) transition, corrected for 
hyperfine structure [31] . 



III. PRECISION MEASUREMENTS IN 
MOLECULES 



High-precision molecular spectroscopy on neutral and 
ionic molecular hydrogen is reviewed in this section. The 
excellent agreement found between these experimental 



results with ab initio calculations, provide the most strin- 
gent tests on the application of quantum electrodynamics 
in a chemically-bound system. 

Recently, the rotationless vibrational transitions {v" = 
—>■ v' = 1) for H2, HD, and D2 were determined to an 
accuracy of ^ 2 x 10^'' cm^^ 44J. These fundamen- 
tal ground tone vibrations of H2, HD, and D2 were ob- 
tained from Doppler-free laser spectroscopy in the colli- 
sionless environment of a molecular beam. The rotation- 
less fundamental vibrational splitting was derived from 
the combination difference between electronic excitation 
from the X^Y.'^^v = and v ^ 1 levels to a common 
EF^Yi'^,v = level. The experimental results are in ex- 
cellent agreement with a full ah initio calculation up to 
an uncertainty 5E ^ 2 x 10^^ cm^^, where the combined 
precision of the experimental and theoretical values is de- 
fined as 5E - 
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Some overtone transition frequencies of molecular hy- 
drogen have been determined via single-photon infrared 
absorptions, with the most recent investigations employ- 
ing cavity-ringdown (CRDS) techniques [41-43 . Despite 
the sensitivity of CRDS, the extremely weak quadrupole 
transitions necessitate high-pressure samples subject to 
pressure shifts in addition to Doppler-broadening. The 
results from Hu et al. [12], Campargue et al. [H] and 
Kassi et al. [13] are in agreement with theory, and re- 
sult in a combined precision of the comparison of 5E ^ 
1 X 10~^ cm~^. Maddaloni et al. [lU^ claim an absolute 
accuracy for the S(0) and S(2) transitions of the funda- 
mental ground tone in D2 at 2 x 10~* cm~^, however, a 
comparison to the theoretical value yields a la discrep- 
ancy. Owing to this large inconsistency from Ref. [iO] . 
bounds derived from these transitions are not included. 

The recent and most accurate experimental determi- 
nations on the dissociation limits of H2 [55] j D2 [37] and 
HD [3S], are based on three energy intervals obtained 
from separate spectroscopic investigations. The ioniza- 
tion potential (IP) is the sum of these energy intervals: 
the first interval is between the ground electronic state 
and the EF d = 0, J state; the second between is the EF 
state to a high-np Rydberg state; and the third is be- 
tween the high-np state to the molecular ion H^ ground 
state X t> = 0. The neutral molecule dissociation limit is 
derived by combining the IP with the accurate theoretical 
calculation of the molecular ion dissociation energy and 
the accurate experimental value of the atomic ionization 
energy. These neutral molecule dissociation limits are in 
excellent agreement with the most accurate ah initio cal- 
culations, again demonstrating the correctness of QED 
evaluations at 5E ^ 1 x 10^'^ cm^^ level. 

High-resolution molecular spectroscopy has been per- 
formed on trapped and cooled HD+ [30, 31j as well. The 
use of sympathetic cooling of HD+ by optically cooled 
Be+ ions reduces the Doppler widths of the molecular 
transitions. In a resonance-enhanced multiphoton disso- 
ciation (REMPD) scheme, the ground state HD+ ion is 
first optically excited to a higher vibrational quantum 



TABLE I: Relevant data from neutral and ionic molecular 
hydrogen transitions used to derive constraints. 
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"5E is the discrepancy between theory and experiment 



state. Thereafter, a second photon further excites the 
ion to a dissociative state leading to a loss of the trapped 
HD+ which can be detected. Excellent agreement be- 
tween experiment and theory is observed for the spectro- 
scopic results on the w" = — > w' = 4 transition [3D] in 
HD+, with the experimental accuracy of the hyperfine- 
less transition at 1.7 x 10~^ cm~^. The theoretical value 
is accurate to 2.3 x 10^^ cni~^ resulting in a combined 



uncertainty of 6E = 1.7 x 10 "^ cm . 

The hyperfine structure is partially-resolved in a recent 
vibrational spectroscopy investigation of the z; == — > 1 
band in HD"*" [31] leading to an improved measurement 
accuracy. However, there is currently a 2a discrep- 
ancy between theory and experiment for the hyperfineless 
transition energy for this fundamental vibrational split- 
ting. It appears plausible that this deviation is caused 
by statistical noise or an as of yet unaccounted hyper- 
fine interaction. Therefore, rather than considering this 
a possible manifestation of a fifth force, we include this 
transition in the fifth force constraining analysis with SE 
set equal to the 2a discrepancy between theory and ex- 
periment. 

Experiments on the w = — )• 8 band in HD+ arc in 
preparation, with great potential for probing fifth force 
interactions, but no precision frequency results have been 
obtained yet [35]. 

The most accurate dataset that is used to further de- 
rive constraints for fifth-force interactions is listed in Ta- 
ble U 



IV. SEARCH FOR FIFTH FORCES IN 
MOLECULES 



The occurrence of fifth forces beyond the SM can be 
phenomenologically parameterized by a Yukawa-type po- 



tential with an effective range A 



(2) 



where /?' is a coupUng strength, which may a priori be 
attractive or repulsive. An observation of a significant ef- 
fect would prompt a more thorough investigation of par- 
ticular model potentials that best represent an observed 
discrepancy. The presence of a Yukawa potential correc- 
tion implies the existence of a force carrier with a mass 
inversely proportional to A. For the effective separation 
distance on the A-scale relevant to this study, this may 
be viewed as searches for effects of new force-carrier par- 
ticles with light masses in the order of keV, where the 
particle interaction length is taken to be the Compton 
wavelength (h/mc) of a hypothetical bosonic gauge par- 
ticle. 

In this particular case of molecules, the effect of a 
fifth-force between nuclei can be searched for at the dis- 
tances where they are bound within the geometry of the 
molecules; in this case r is the internuclear distance. To 
incorporate the effects of different nucleon numbers N 
in H2, HD, and D2, as well as the Hj and HD+ ions 
a redefinition of the coupling constant is introduced to 
explicitly express the dependence: 



V5{r)=l3NiN2Y{r) 



(3) 



where A^i_2 are the nucleon numbers for each nucleus; 
note that this differs from a definition by Bordag et 
al. |53j . The extra long-range hadron-hadron interac- 
tion probed here is spin- independent. In view of the nu- 
cleon scaling factors, molecules with the highest number 
of nucleons will provide the tightest constraint on the 
existence of fifth forces parameterized by V5(r); hence 
the D2 isotopomer would be a more sensitive test ground 
than H2. 

We treat the extra potential in Eq. pi as a perturba- 
tion on the ro- vibrational level energies of the molecular 
states. For a transition between the ground v" and ex- 
cited v' vibrational levels, the contribution (AV5) of a 
fifth-force potential can be expressed as 

(AV^s.a) = /?iViiV2[(*.-,j'(r)|y(r,A)|*,,,j,(r)) 

-(*.",j"W|i"(^,A)|*,»,j"W>] 
= /37ViiV2An (4) 

where "^y^j^r) represents the nuclear wave functions. Al- 
though included in the calculations, the explicit mention 
of rotational quantum numbers J is omitted in the rest of 
the discussion since rotation has a much smaller contri- 
bution compared to vibration for the transitions treated 
here. 

Numerical calculations were performed to evaluate the 
(AVs^a) contribution of various electronic ground state 
vibrational transitions using the accurate wave functions 
for H2, HD, and D2 [32H31] and HD+ [33. The differ- 
ence of the expectation values AYx in Eq. (|4|) is plotted in 
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FIG. 1: (Color online) Calculated difference of the ex- 
pectation values AY\ for different values of A for several 
{v' ,v" — 0) vibrational energy separations in H2, D2, HD 
and HD"*" molecules. 



Fig. [T]for the w" = — > f' = 1, 2 transitions in H2, with 
the interaction length A taken as a parameter. The cor- 
responding differential contribution for vibrational tran- 
sitions for v" — ^ v' = 1 for D2 and HD and for 
w" = -^ u' = 1, 4 in HD+ are also plotted in Fig. [ij 

This illustrates that for a specific molecule, the sensi- 
tivity for probing a fifth-force contribution increases as 
Av increases. This can be attributed to the difference 
in the spatial extent of the ground and excited state 
wave functions, where the latter has a most probable 
position (related to classical bond length) gradually dis- 
placed from the equilibrium for increasing v quantum 
numbers. Since the wave function density for consecu- 
tive vibrational levels shifts only slightly, of effect of a 
Yukawa potential VJ5 is probed in molecules as a differen- 
tial contribution. For a given vibrational transition, H2 
has the greater sensitivity compared to D2 and HD (be- 
fore the nucleon-number scaling) since the wave functions 
belonging to different vibrational levels have more simi- 
lar spatial extents in the heavier isotopes than in H2 . For 
specific vibrational transitions, e.g. {v" = ^- w' = 1), 
the differential contribution for H2 is also greater than 
the corresponding contribution HD+ since the distance 
between the nuclei in the ions, with only one electron 
contributing to the chemical bond, is larger. In general, 
transitions in molecules with shorter bond lengths will be 
more sensitive to a fifth-force hadron-hadron potential. 

Any difference SE between the experimental and cal- 
culated values for a particular transition energy can be 
used to set bounds for the maximum contribution to a 
fifth-force. A constraint for the coupling strength f3 is 
obtained for a range of values of an interaction length A 
by the relation 



/?< 
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N1N2AYX ' 



(5) 



The coupling constant /3 has units of energy x distance. 



and can also be related to any known interaction such 
as the electromagnetic interaction characterized by the 
coupling constant a. The Coulomb potential can be ex- 
pressed as Vem = otZxZ^jr with units of energy/distance, 
where Z\^i — 1 for the nuclei considered in this study. 
Except for the extra exponential factor, V5 has the same 
1/r-form and hence the same dimensions as Vem, lead- 
ing to a dimensionless ratio of their interaction strengths 
/3/a. 



The fundamental vibration in the hydrogen 
molecule 



The agreement of the accurate experimental and the- 
oretical values are used to provide a constraint on the 
effect of new interactions. The combined experimen- 
tal and theoretical uncertainty at &E = 2 x 10^'* cm^^ 
for the fundamental vibrational tone provides limits on 
the interaction strength |/3| for different values of the 
interaction length A. The limits derived from the H2 
X, w" = — >■ ri' = 1 transition constrain the strength 
of a new interaction to be /3 < 4.7 x lO^^a for inter- 
action lengths A > 1 A. Similarly bounds are obtained 
using the HD and D2 fundamental vibrational tone at 
/3 < 3.4 X lO^^a and /3 < 1.5 x 10~*a, respectively, for 
interaction lengths A > 1 A. 



B. Overtone vibrations in neutral hydrogen 

Similar constraints can be derived from investigations 
of the electronic ground state direct overtone vibrations 
(e.g. w' = to w = 2,3) of H2 and D2. The results 
from recent CRDS studies on the overtone quadrupole 
transitions of H2 (HI gS] and D 2 [101 US] can be used 
to extract constraints. The overtone transitions of H2 
and D2 have intrinsically higher sensitivity of the tran- 
sitions in comparison to the fundamental ground tone as 
shown in Fig. 2] However, the worse uncertainty (^ 10~'^ 
cm~^) from these Doppler-limited studies, and pressure- 
shift corrected transition energies, result in less tight con- 
straints. For interaction lengths A > 1 A, the H2 (2,0) 
band constrains /3 < 4.7 x 10~^a and the (3,0) band con- 
strains /3 < 2.4 X 10~^a, while the D2 (2,0) band leads 
to a constraint of /3 < 4.2 x lO^^a. 



C. Level energies in the HD^ molecular ion 

The spectroscopic results for HD+ for the w" = -> 
ti' = 4 transition |30j provide stringent bounds on possi- 
ble fifth-force interactions. For A > 1 A, the constraint 
for the interaction strength is /3 < 1.1 x 10~^a. This 
tighter bound from HD"*" is due to the better accuracy in 
both experiment and theory, as well as larger Aw probed 
between the w = and w = 4 quantum states. 
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FIG. 2: (Color online) Calculated expectation value Yx^dq of 
a fifth-force contribution to the H2 dissociation limit Do for 
a range of A values. For reference the contributions AIa for 
the {v' = 1, v" = 0) band of H2 and the (y' = 8, v" = 0) band 
of HD"*" electronic ground state are plotted as well. 



The results for HD+ vibrational spectroscopy of the 
z; = — > 1 transitions can still be used to also pro- 
vide bounds for new hadronic interactions. The two- 
sigma discrepancy between theory and experiment 8E = 
4.7 X 10^^ cm^^ is used to derive bounds to a fifth-force 
interaction with /3 < 1.1 x lO^^a for A > 1 A. 



D. The dissociation limit in H2, HD and D2 

The dissociation energy is defined to be the energy dif- 
ference between the deepest bound molecular state v" = 
0, J" = and the state when the two constituent atoms 
are non-interacting, i.e. ai r — 00 where V(oo) = 0. In 
a similar way, a possible fifth-force interaction is treated 
perturbatively and the contribution is expressed as 



(AV. 



5, a; 






.o(r)|r(r,A)|M'.».o(r)) 



(6) 



where A is again treated as a parameter in the calcula- 
tions. 

The fifth force contribution to the dissociation limit is 
the upper limit of the Aw-progression for the vibrational 
transitions represented in Eq. Q, and thus intrinsically 
is more sensitive than any of the ground state vibrational 
splitting. 

The combined experiment-theory accuracy of Dq for 
H2 is 5E = 1.2 X 10-3 cm-i Ref. [35]; for D2 [37] it is 
5E = 1.1 X 10^3 cm-i; and for HD [35] it is 5E = 1.1 x 
lO^"^ cm^^. The D2 value gives the tightest constraint 
at /3 < 3.8 X 10~^a for interaction lengths A > 1 A, H2 
and HD gives ;3 < 1.7 x lO^^a and /3 < 7.5 x lO^^a, 
respectively, for A > 1 A. 
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FIG. 3: (Color online) Limits on the strength of the coupling 
constant /3 relative to the electromagnetic coupling constant 
Of as a function of the interaction range A. The alternate axis 
expresses /3 in units of eV-A. 



The present bounds obtained from the D2 dissocia- 
tion limit is only a factor four less tight compared to 
that of the HD+ (4,0) band, despite the experimental 
uncertainty of D2 Dq being ^60 times worse than for the 
ion. This suggests that future improvements in the ac- 
curacy of the D2 Dq [54] can lead to tighter constraints, 
while in the case of HD+ hyperfine interaction need to 
be addressed more accurately to improve its accuracy. In 
addition, as the wave functions of the neutral hydrogen 
molecules extend towards shorter range than the ionic 
species, the bounds for /3 at shorter A is tighter for the 
neutrals, e.g. for A < 0.4 A in Fig. |3|l?o of D2 gives the 
tightest constraints. 

Precision measurements on exotic molecules, such as 
anti-protonic helium |55j . could provide similar con- 
straints for force ranges on shorter length scales, i.e. 
A < 0.5 A. 



VI. CONCLUSION 



V. DISCUSSION 

Constraints on a fifth-force coupling constant /3, for a 
phenomenological Yukawa interaction potential Vc,(r) — 
/3exp(— r/A)/r, are derived from vibrational transitions 
in both neutral and ionic molecular hydrogen species 
measured at high precision. Some of the tightest con- 
straints are plotted in Fig.lsj where the value of \(3\/a is 
used to express the strength of a fifth-force with respect 
to the strength a of the electromagnetic interaction, with 
an alternative axis expressing /3 in units of eV-A. The re- 
gion below a curve represents the allowed value of /? for 
a certain value of the range of the force A. The con- 
straints obtained from HD+ are most stringent owing to 
the better accuracy in both the experimental and calcu- 
lated values. At present, the HD+ (4,0) band furnishes 
slightly tighter limits than that of the (1,0) band, both 
because of the inherent enhancement {AY\) of the for- 
mer and also the 2a deviation between theory and ex- 
periment in the latterwhich might be due to a statistical 
error or to an as of yet unaccounted hyperfine interac- 
tion. The tighter constraint from D2 compared to H2, 
gained from the nucleon number scaling, is demonstrated 
graphically in Fig. [31 The present analysis yields a con- 
sistent constraint P/a < 10~® or /3 < 1 x 10~^ eV-A for 
long-range hadron-hadron interactions at A > 1 A. The 
allowed value for /3 is indicated by the shaded region in 
Fig.[3| 



The advancements in recent years on the ab initio the- 
ory of light molecular species, e.g. neutral and ionic 
molecular hydrogen and its isotopomers, have led to 
the successful application of QED corrections in quan- 
tum chemical calculations. Comparison of these theo- 
retical results with highly accurate experimental values 
from precision spectroscopic investigations have demon- 
strated excellent agreement between theory and experi- 
ment. Treating fifth force long-range hadron-hadron in- 
teractions in the form of a Yukawa-type interaction po- 
tential, we show how such highly accurate comparisons 
between theory and experiment provide a search ground 
for new physics beyond the Standard Model. 
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